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Abstract 

We measure Bose-Einstein correlations between like-sign charged pion pairs in ha- 
dronic Z decays with the L3 detector at lep. The analysis is performed in three 
dimensions in the longitudinal center-of-mass system. The pion source is found to 
be elongated along the thrust axis with a ratio of transverse to longitudinal radius 
of 0.81 ± 0.02+°;?3. 
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1 Introduction 



It has long been realized that the shape and size in space-time of a source of pions can be 
determined, as a consequence of the interference of identical bosons, from the shape and size of 
the correlation function of two identical pions in energy- momentum space The space-time 

shape of a source during hadronization is important experimental information on QCD in a 
sector where perturbative methods are not applicable. 

The form of the correlation function in more than one dimension has been a major subject 
of theoretical study in recent years [p|-[TO[]. In Monte Carlo generators, spherical symmetry is 
usually assumed |[TT|-[15|, while elongation can be expected when a string-like shape is main- 
tained [p|,[lO,|16|. Experimentally, detailed three-dimensional analyses have been done only for 



heavy ion collisions |17,|T^ and for hadron-hadron collisions |19|. While the volume of the pion 
emission region appears to be approximately spherical for heavy ion collisions, a clear elonga- 
tion is observed in hadron-hadron collisions. Analyses in e"'"e~ collisions have generally been 
limited to one dimension |^. An exception is a preliminary analysis |^|, which indicates an 
elongation in e"'"e~ collisions. 

Recently there has been a revival of interest in Bose-Einstein (BE) correlations of particles 
produced in the hadronization of energetic quarks, mainly due to its possible impact on the 
measurements of the W mass UTT Hl^] in the four-jet channel e"'"e~— >W+W~— > qqqq. The ex- 
perimentally observed radius of hadron emission is of order 1 fm, about an order of magnitude 
larger than the distance between the W''"W~ decay vertices at present energies. The actual 
overlap in W"'"W~— > qqqq decays is of course determined by the (possibly non-spherical) con- 
figuration of the decay partons of the W's in all space-time dimensions during hadronization. 
A better understanding of BE correlations in e+e" — Z — qq, using high statistics lep data, 
will contribute to understanding the BE effect in W'^W". 

In this letter, two-particle correlations of like-charged pions are studied in hadronic Z-decay 
as a function of three components of the four-momentum difference Q. Effects due to hard 
gluon radiation, Z-decay into heavy quarks or the dependence on the transverse mass of the 
particles are not considered. 



2 Data 

2.1 Charged-Particle and Event Selection 



The data used in the analysis were collected by the L3 detector in 1994 at a center-of-mass 



energy of y/s ~ 91.2 GeV. The data selection uses information on charged particles from the 



Time Expansion Chamber (tec) and the Silicon Microvertex Detector (smd) |22]. 

To obtain a data sample of hadronic Z decays, we perform an event selection using charged 
tracks. The charged tracks are required to have at least 40 (of 62 possible) hits in tec, and 
the number of wires between the first and last hit is required to be at least 50. The distance 
of closest approach (projected onto the transverse plane) of a track to the nominal interaction 
vertex is required to be less than 5 mm. The transverse momentum of a track must be greater 
than 100 MeV. 

In order to reduce background arising from beam-gas and beam-wall interactions as well as 
from leptonic events and from two-photon interactions, we use the following criteria 

Edfi>0.15, Ei^<0.75, S^<o.75, Ar,h>4, (1) 
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where Pi is the momentum of particle i, with components p\\i parallel to the beam direction and 
in the plane perpendicular to the beam direction, and where the sums run over all the tracks 
of the event. N^h is the number of charged-particle tracks. A small remaining background of 
T^T~ events in which both r particles decay into more than one charged particle is removed by 
requiring the second largest angle (^2 between any two neighboring tracks in the R — ip plane 
to be less than 170°. After this cut the r+r^ contamination is approximately 8 ■ 10~^. 

To make sure that events lie within the full acceptance of the tec and SMD, we require 
|cos^thr| < 0.7, where ^thr is the polar angle of the thrust axis determined from charged tracks 
only. 

A total of about one million events satisfy the selection criteria. 

For the computation of the four-momentum difference Q, the resolution of the angle between 
pairs of tracks is of crucial importance, especially for small Q- values. For this reason we impose 
an additional cut, requiring an unambiguous polar angle measurement. This ensures good 
resolution of variables, such as Q, which depend on two tracks. Together with the previous 
cuts about 40% of the tracks are rejected. The final data sample consists of approximately 36 
million like-sign track pairs. 

With this selection, good agreement is obtained between data and simulation for the distri- 
butions of the differences between pairs of tracks of the azimuthal and polar angles with respect 
to the beam and for the distributions of and the components of Q used in this analysis. 
This is shown in Fig. ||, where the data are compared to events generated with jetset 



(including BE effects) which have been passed through the L3 detector simulation program |2i 
reconstructed, and subjected to track and event selection in the same way as the data. Simi- 



lar comparisons using events generated by jetset without BE effects or by herwig also 
show good agreement except at small values of the variables, where the Bose-Einstein effect is 
important. 



3 Analysis 

3.1 Bose-Einstein Correlation Function 

The two-particle correlation function of two particles with four-momenta pi and p2 is given by 
the ratio of the two-particle number density, P2{Pi,P2), to the product of the two single-particle 
number densities, Pi{pi) Pi{p2) ■ Since we are here interested only in the correlation R2 due to 
Bose-Einstein interference, the product of single-particle densities is replaced by ^0(^1)^2)5 the 
two-particle density that would occur in the absence of Bose-Einstein correlations: 

R2{Pi,P2) = —. r • (2) 

Po{PuP2) 

i?2 — 1 is related to the space-time particle density through a Fourier transform 

Since the mass of the two identical particles of the pair is fixed to the pion mass, the 
correlation function is defined in six- dimensional momentum space. Since Bose-Einstein cor- 



relations can be large only at small four-momentum difference Q = ^/—{pl — ^2)^; they are 
often parametrized in this one- dimensional distance measure. There is no reason, however, to 
expect the hadron source to be spherically symmetric in jet fragmentation. This is the reason 
for performing a three-dimensional analysis. 
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3.2 Longitudinal Center-of-Mass System 

In our analysis we use the longitudinal center-of-mass system (LCMS) f^. This is defined for 
each pair of particles as the system, resulting from a boost along the thrust axis, in which the 
sum of the momenta of the pair is perpendicular to the thrust axis. In this system, we can 
resolve the three-momentum difference of the pair of particles into a longitudinal component 
Ql parallel to the thrust axis, Qout along the sum of the particles' momenta (see Fig. 0) and 
Qside perpendicular to both Ql and Qout- Then, the invariant four- momentum difference can 
be written as ^ 

Q' = Ql + Qsidc + Qlut - im' = Ql + Qsidc + QL(i - P'), (3) 

where 

n _ Pout 1 + Pout 2 / , \ 

with Pouti and Ei {i = 1, 2) the out-component of the momentum and the energy, respectively, 
of particle i in the LCMS. The energy difference AE and therefore the difference in emission 
time of the two particles couples only to the component Qout- Consequently, Ql and Qside reflect 
only spatial dimensions of the source, whereas Qout reflects a mixture of spatial and temporal 
dimensions. The correlation function is then parametrized in terms of Q = {Ql, Qside, Qout)'- 

= . (5) 

poiQ) 

3.3 Determination of the Correlation Function R2 

After determining p2{Q), there are three steps in obtaining R2. The first step is to create a 
'reference sample', which is done by event mixing, and from it to determine Po{Q)- Then two 
corrections must be determined, one for non-BE correlations lost in the event mixing and one 
for detector effects. 

The reference sample, from which po is determined, is formed by mixing particles from 
different data events in the following way. First, events are rotated to a system with the 2;-axis 
along the thrust axis and are stored in a 'pool'. Then, events are randomly selected from this 
pool and their tracks replaced by tracks of the same charge from other events in the pool of 
approximately the same multiplicity under the condition that no track originates from the same 
event. Finally, Ql, Qout and Qside are calculated for each pair. Used events are removed from 
the pool, thus preventing any regularities in the reference sample. From this mixed sample we 
obtain the particle density Pmix(Q)- 

Since this mixing procedure removes correlations other than just those of Bose-Einstein, 
e.g., those from energy- momentum conservation and from resonances, pmix is then corrected 
for this loss by a factor Cmix, which is estimated by Monte Carlo (MC) using a generator with 
no Bose-Einstein (BE) effects (jetset or herwig). Thus, in the absence of Bose-Einstein 
correlations, the corrected two-particle density is given by 



PoiQ) = Pmix(<5) ■ C'mix(Q) , where C^miQ) 



P2{Q) 
Pmix(<5) 



(6) 

MC.noBE 



The ratio (p2/pmix) must further be corrected for detector resolution, acceptance and effi- 
ciency and for particle misidentification. For this we use a multiplicative factor derived from 
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Monte Carlo studies. Since the L3 detector does not identify the hadrons, this factor, C^et, is 
given by the ratio of the two-pion correlation function found from MC events at generator level 
to the two-particle correlation function found using all particles after full detector simulation, 
reconstruction and selection: 



Cdet(Q) 



P2{Q) 
Pmix('5) 



gon, pions 




P2iQ) \ 

x(Q)/ 



(7) 



det, all 



Taking all charged particles, instead of only pions, in the generator level MC, leads to consistent 
results. Combining this correction factor with (^) and results in 



R2iQ) 



P2{Q) 1 

Pmix (Q) c^ix(g) 



(8) 



The analysis is done in three-dimensional bins of Q. In terms of numbers Nkim of like-sign 
particle pairs in the three-dimensional bin k,l,m of Ql, Qont and Qside, eq. (H) becomes 



R 



■2 klm 



klm 



^^klm 



data 



Armix 
^^klm 



klm 



MC,noBE 



^^klm 
Ti rgen, mix 
^^klm 



n: 



det , mix 
klm 



jvrdet 
^^klm 



(9) 



MC 



Each Nkim is normalized to the total number of pairs in the corresponding sample. 

The resolution in the variables Qi is estimated using Monte Carlo events to be 0.02-0.05 
GeV for Qi < 0.2 GeV. Given the available statistics, we choose a bin size of 0.08 GeV. 

In our analysis, we use jetset without BE and herwig (not having BE) to determine 
the mixing correction factor Cmix and jetset (with and without BE) as well as herwig to 
determine the detector correction factor Cdct- Together with a variation of the mixing technique, 
the selection criteria and the fit range, these six Monte Carlo combinations will serve to estimate 
systematic errors. 



3.4 Parametrization of R2 

Assuming a Gaussian (azimuthally, but not necessarily spherically, symmetric) shape of the 
source, the following three-dimensional parametrization has been proposed [p|,|6|,p6|: 

R2{Ql, Q out, Q side) = 

7 (1 + 5Ql + eQont + ^Qside) 

1 + X exp [-RIQI - -RoutQout - -^sidcQsidc + 2pL,out-RL-Rout<5LQout) , (10) 

where the factor (1 + 6Q1, + eQont + ^<5sido) takes into account possible long-range momentum 
correlations in the form of a slow rise, 7 is a normalization factor close to unity and the term 
between square brackets is the two-particle Bose-Einstein correlation function associated with 
a Gaussian shape of the source. 

By fitting the correlation function with this parametrization, one can extract the incoherence 
factor A, which measures the strength of the correlation, and the 'radii' Ri {i = L, out and 
side) defined as l/\^ai, with af the variances of a multi-dimensional Gaussian distribution 
of the source in configuration space. pL,out is the correlation between the longitudinal and out 
components of this Gaussian. In the LCMS, the duration of particle emission only couples to the 
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out-direction and only enters in the parameters i?out and pL.out- Hence, -Rside can be interpreted 
as the transverse component of the geometric radius. The parametrization, eq. (|l3), assumes 
azimuthal symmetry of the source, which means that the two-particle Bose-Einstein correlation 
function associated with the Gaussian shape of the source, is invariant under the transformation 
Qside —Qside- Consequently, the only possible off-diagonal term is the QL^out term. 

We first checked the method using Monte Carlo events at detector level, which were gener- 
ated by JETSET without BE, instead of data. The correction factor Cmm was determined using 
JETSET without BE and using herwig, while Cdet was determined using jetset with BE and 
using HERWIG, giving four different combinations of correction factors. Also we used events 
generated by herwig as data, with the same two MC determinations of Cmix but with jetset 
with BE and jetset without BE for Cdet- In all eight cases fits gave results consistent with 
A = (x^/NDF ~ 1 for A = 0), as expected in the absence of Bose-Einstein correlations. 

4 Results 

The results of a three-dimensional fit, in the range Qi < 1.04 GeV, of the Bose-Einstein correla- 
tion function R2 with the parametrization of eq. (p!oD are presented in Table 0. The off-diagonal 
term turns out to be zero within errors (pL.out = —0.008 ± 0.057) and the results given in the 
table correspond to a fit with this term fixed to zero. The values are obtained using jetset 
without BE for the mixing correction and JETSET with BE for the detector correction, since 
the latter model is found to be the most successful in reproducing our data in the relevant 
variables (Fig. p and since this choice gives the lowest in the fits. 

To estimate the systematic errors on the fit parameters, we examined four different sources, 
repeating the full analysis in each case. First, we looked at the fit results obtained with the six 
possible combinations of mixing and detector Monte Carlo corrections. The systematic error 
from this source is taken as the RMS of these six values. Secondly, the influence of a different 
mixing sample was studied by removing the conditions that tracks are replaced by tracks with 
the same charge and coming from events with approximately the same multiplicity. For each of 
the six Monte Carlo combinations the difference in the results between the two mixing methods 
was taken as an estimate of the systematic error from this source and the square root of the 
mean of the squares of these differences taken as the systematic error from this source. In the 
same way systematic errors related to track/event selection and to the choice of fit range were 
evaluated. The analysis was repeated with stronger and weaker selection criteria, resulting in 
approximately 11% fewer/more events and 6% fewer/more tracks. The upper limit of the fit 
range was varied by ±0.16 GeV in all Qi. As a total systematic error (second error in the table) 
we add the four errors from these sources in quadrature. 

For i?L the errors from all the sources are approximately equal. However, the systematic 
error on A and the transverse radii is dominated by the contribution from the six possible MC 
combinations to determine Cdot • C'mix- Of the two, Cdet shows the larger generator dependence 
and accounts for most of the systematic error. It is approximately unity for most values of Q, 
but increases to about 1.1 near Q = in the case of jetset with BE and decreases to about 
0.95 for JETSET without BE and to about 0.9 for herwig. We understand this difference in 
behavior between jetset with BE and the two Monte Carlos without BE as follows. Since Q 
depends both on the energies of the particles as well as on the angle between them, small Q can 
be due to either small angle or low energies. In a Monte Carlo with BE the fraction of pairs 
at small Q with small angle will be larger than in the other Monte Carlo models. This will 
lead to lower detection efficiency and hence to larger corrections. This part of the error is also 
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responsible for the asymmetry in the errors since our best choice of MC combination does not 
coincide with the average. In particular, all of the other combinations result in smaller values 
of -Rside/ Rl- 

To study the behavior of the components of Q, projections of R2, eq. @, onto the three 
axes are shown in Fig. using the regions Qi < 240 MeV {i.e., the first three bins) of the 
non-projected components. The dashed curves correspond to the fit results described above. 
Similarly, a projection of R2 onto the Qside-^L plane is shown in Fig. ^ From Table |I| we find 
that i?L is larger than both i?out and -Rside- The ratios Rside/ Rl and Rout/Rh are 5 standard 
deviations less than unity. Thus, the pion source is elongated along the thrust axis. 

From the value of (see Table ||) it appears that the shape of the correlation function 
deviates from a Gaussian. An exponential has often been suggested as an alternative to the 
Gaussian. Accordingly, the fits were repeated with an exponential, exp(— i^jQj), replacing 
the Gaussian in eq. ( |T0|) with pL.out = (dotted lines in Fig. |]). Although the overall 
increases by typically 2-4%, depending on the correction combination used, the fits confirm the 
elongation observed from the Gaussian fit. 

A general approach to study deviations from the Gaussian is to use an expansion, 
due to Edgeworth p8[, in terms of derivatives of the Gaussian, which are related to Hermite 
polynomials. Taking only the lowest-order non-Gaussian term into account, results in 

-R2(Ql) Qout) Qsidc) = 

7 (1 -h + eQont + ^Qside) 

■ |l + Aexp [-RIQI - -RoutQout - RlideQlide) 

■ 1 + -^Hsi^Ri^Qi,) ■ 1 + -^Hsi^RontQont) ■ 1 H ^-f^3(-RsidcQside) | , (H) 

where Ki {i = L, out, side) is the third-order cumulant moment in the corresponding direction 
and H^{RiQi) = {\/2RiQiY — 3\/2RiQi is the third-order Hermite polynomial. Note that the 
second-order cumulant corresponds to the radius Ri. The results of the fit with the Edgeworth 
expansion are shown in Table |l]and Fig. ^ (full lines). The value of x^/NDF and the confidence 
level indicate a better fit than the Gaussian one. The non-zero values of the k parameters 
indicate the deviation from a Gaussian. A is larger than the corresponding Gaussian A. The 
values of the radii confirm the elongation observed from the Gaussian fit. 

For comparison we have also performed fits on a JETSET generator sample taking R2 as 

p _ (A^/A^™"") JETSET with BE 

''''^^''^^ " (Ar/Ar--)jETSETnoBE ' 

The resulting values of Rside/Rh and Ront/Rh are 1.08±0.03 and 0.79±0.02, respectively. This 
value of Rside/Rh is larger than in the data (see Table |1|). Thus, the standard Bose- Einstein 
implementation of jetsetQ fails to reproduce the experimentally observed elongation. It is 
worth noting that even though there is no explicit spatial asymmetry in the JETSET treatment 
of BE, neither -Rside nor -Rout is found to be equal to -Rl. 

Another Monte Carlo generator with Bose-Einstein simulation is VNI |^, which in its 
present form has been found to predict -Rside ~ -Rl [0) in contradiction to our results. 



Since the difference between the two transverse components of -R is small compared to that 
between the longitudinal component and either of the transverse components, we check the 



*^The Bose-Einstein simulation is done by the subroutine luboei, with the l3 defauh values parj(92)=1.5 
and PARj(93)=0.33 GeV. 
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results of the three-dimensional analysis by a two-dimensional analysis in i?L and -Rt, where we 
can use intervals of 40 instead of 80 MeV. The out and side terms in the exponential of eq. ([T0|) 
are replaced by R'^Q'^, with = Qlnt + Qside- The two-dimensional fits (both Gaussian and 
Edgeworth) result in values of A consistent with those from the three-dimensional fits and with 
values of Rt/Ri, lying between the values of Rout/Rh and Rside/ Rl, confirming the elongation 
observed in the three-dimensional fits. As in the three-dimensional fits, the confidence level 
of the Gaussian fit is poor, that of the exponential fit poorer, and that of the Edgeworth fit 
acceptable. These statements are also true for analyses of jetset. 



5 Summary 

A sample of 1 million events of the electron-positron annihilation process e"'"e~ — Z — » hadrons 
at LEP, has been used to study two-particle Bose- Einstein correlations of like-charged pions. 
The analysis was performed in three dimensions, as well as in two dimensions, in the longitu- 
dinal center-of-mass system. Using the Gaussian parametrization of eq. (ITUp we extracted the 



incoherence factor. A, and the radii, i?L, -Rout and -Rside- A better fit was obtained using the 
Edgeworth expansion (|TT|) of the Gaussian parametrization. The transverse radius is found to 
be significantly smaller than the longitudinal radius: 



Rside/ Rl = 0.81 ± o.o2i°: 



0.03 

19 



Our measurement implies that models based on the assumption of a spherical source are 
too simple. 
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Table 1: Values of the fit parameters for the three-dimensional analyses of L3 data. The first 
error is statistical, the second systematic. 



parameter 


Gaussian 


Edgeworth 


A 


0.41 ± OMf^il 


0.54 ± 0.02lg:^^ 


Rl (fm) 


0.74 ± 0.02to:o3 


0.69 ± 0.02lg:°^ 


-Rout(fm) 


0.53 ± 0.021^:^^ 


0.44 ± 0.02lg:g^ 


-Rside (fm) 


0.59 ± OMtoil 


0.56 ± 0.021°:°^ 




0.71 ± 0.02+°:°! 


0.65 ± 0.03lg:g^ 


-Rside/ Rl 


0.80 ± 0.02lg:?| 


0.81 ± 0.02lg:?^ 


«L 




o.5±o.ii°i 






0.8 ±0.1 ±0.3 


'^side 




0.1 ±0.1 ±0.3 


6 


0.025 ± 0.005ig:gJ^ 


0.036 ± omitoml 


e 


0.005 ± 0.005lg:g?2 


0.011 ± 0.005^^:^?^ 




-0.035 ± 0.005lg:g^^ 


-0.022 ± 0.006lg:g^^ 


xVndf 


2314/2189 


2220/2186 


C.L. (%) 


3.1 


30 
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Figure 1: Distributions of (a) the difference in polar angle of pairs of tracks, 59^ (b) the difference 
in azimuthal angle of pairs of tracks, 50, (c) the four- momentum difference squared, Q^, and the 
(d) longitudinal, (e) out, and (f) side components of the four-momentum difference Q, for data 
(points) compared to the predictions of jetset with BE after detector simulation (histogram). 
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Pi 




Figure 2: The longitudinal center of mass frame (LCMS) showing the projection of Q on the 
{Qh-Qout) plane. Qside is the projection of Q on the axis perpendicular to this plane. 
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Figure 3: Projections of i?2 onto the three axes Ql, Qout and Qside using the regions up to 
240 MeV of the non-projected components. The mixing correction is determined using jetset 
without BE and the detector correction using jetset with BE. The full lines correspond to 
projections of the fit with the lowest-order Edgeworth expansion, the dashed lines to those of 
the Gaussian fit and the dotted lines to the exponential fit. 
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Figure 4: Projection of R2 onto the Qside-<5L plane using the regions up to 240 MeV of Qout- 
The mixing correction is determined using JETSET without BE and the detector correction 
using JETSET with BE. 
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